† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11304079, 11304140, 11404094, and 11504088), the China National Scholarship Foundation (Grant No. 201508410255), the Foundation for Young Core Teachers of Higher Education Institutions of Henan Province of China, the Foundation for Young Core Teachers of Henan University of Technology in China, the Korea Institute of Science and Technology (KIST) Institutional Program (Grant No. 2E26291) and Flag Program (Grant No. 2E26300), and the Research Grants of NRF funded by the National Research Foundation under the Ministry of Science, ICT & Future, Korea (Grant No. NRF-2015H1D3A1036078).
Hydrogen spillover mechanism of metal-supported covalent-organic frameworks COF-105 is investigated by means of the density functional theory, and the effects of metal catalysts M4 (Pt4, Pd4, and Ni4) on the whole spillover process are systematically analyzed. These three metal catalysts exhibit several similar phenomena: (i) they prefer to deposit on the tetra (4-dihydroxyborylphenyl) silane (TBPS) cluster with surface-contacted configuration; (ii) only the H atoms at the bridge site can migrate to 2,3,6,7,10,11-hexahydroxy triphenylene (HHTP) and TBPS surfaces, and the migration process is an endothermic reaction and not stable; (iii) the introduction of M4 catalyst can greatly reduce the diffusion energy barrier of H atoms, which makes it easier for the H atoms to diffuse on the substrate surface. Differently, all of the H2 molecules spontaneously dissociate into H atoms onto Pt4 and Pd4 clusters. However, the adsorbed H2 molecules on Ni4 cluster show two types of adsorption states: one activated state with stretched H–H bond length of 0.88 Å via the Kubas interaction and five dissociated states with separated hydrogen atoms. Among all the M4 catalysts, the orders of the binding energy of M4 deposited on the substrate and average chemisorption energy per H2 molecule are Pt4 >Ni4 >Pd4. On the contrary, the orders of the migration and diffusion barriers of H atoms are Pt4 <Ni4 <Pd4, which indicates that Pt4 is the most promising catalyst for the hydrogen spillover with the lowest migration and diffusion energy barriers. However, the migration of H atoms from Pt4 toward the substrate is still endothermic. Thus direct migration of H atom from metal catalyst toward the substrate is thermodynamically unfavorable.
Demand for energy sources is dramatically increasing with the rapid increase of population, while the traditional fossil energy is limited and the use of fossil energy leads to the emission of greenhouse gases.[1] This conflict has led to the increasingly urgent need for clean and ecologically friendly sources of energy.[2] Hydrogen has been regarded as a promising candidate due to its inherent advantages such as high energy density, clean, safe, sustainable, etc.[3] However, hydrogen storage remains one of the main challenges to the realization of fuel cell powered vehicles using hydrogen as an energy carrier.[4,5] Therefore, many different approaches and technologies for hydrogen storage have emerged, in which hydrogen storage on porous materials is one contender. Among various porous materials, covalent-organic framework (COF) is a class of highly porous material, because of its large surface area, extremely low density and porosity, hydrogen storage in COF has attracted significant attention since it was first synthesized in 2005.[6] It has been demonstrated that COF shows promising performances of hydrogen storage at very low temperatures.[7–9] However, its hydrogen uptake dramatically decreases at ambient temperature, due to the very weak van der Waals interactions between H2 and COF (about 2 kJ/mol–5 kJ/mol),[10] which cannot satisfy the requirements for practical applications. It is worth remembering that the binding energy for the hydrogen storage system at moderate conditions (−30 °C–80 °C and 1 bar–100 bar, 1 bar = 105 Pa) should be 9.6 kJ/mol–58 kJ/mol to meet the requirements for practical applications.[11]
Hydrogen spillover has been experimentally demonstrated as an effective approach to improve the hydrogen storage properties of porous materials at room temperature, including activated carbons,[12] carbon nanotubes,[13] zeolites,[14] graphene,[15] and MOFs.[16] Hydrogen spillover is defined as the dissociative chemisorption of H2 molecules on metal catalysts, the H migration onto adjacent surfaces of a receptor and the H diffusion on the substrate surface. Currently, hydrogen spillover has also been successfully applied to COFs.[17,18] For example, Li and Yang[17] studied the hydrogen storage of COF-1 by building carbon bridges, and found that the hydrogen uptakes were enhanced significantly by a factor of 2.6 at 298 K and 10 MPa. Kalidindi et al.[18] showed that the H2 capacities were enhanced by a factor of 2–3 through Pd impregnation on COF-102 at room temperature and 20 bar.
Despite the fact that these promising results have been achieved, the plausibility of spillover mechanism and its usefulness as a hydrogen storage mechanism is still controversial. These controversies mainly arise from the limited understanding of the process that can be afforded through experiment.[19] Theoretical computational work can be very useful in this case. Therefore, many theoretical investigations have also been done on hydrogen spillover in COFs.[8,9,20–22] For example, Suri et al.[20] and Ganz and Dornfield[21] predicted the saturation hydrogen storage capacities of COF-1 and COF-5 by constructing the saturation models. However, the microscopic mechanism of improving the hydrogen storage performance was not involved. Ganz and Dornfield[22] further studied the thermodynamics and energetics of the initial stages of hydrogen spillover on the COF-5 material. They pointed out that spillover could not proceed on pure nor doped COF-5 because of high kinetic barrier. Our group[8] investigated molecular hydrogen and spillover hydrogen storage on several two-dimensional (2D) COFs, and found that their saturation densities were much larger than the physisorption uptakes of H2 molecules at room temperature. Guo et al.[9] used the density functional theory (DFT) to study the hydrogen spillover mechanism of Pt4 supported COFs. They pointed out that although the hydrogen storage capacity of COFs can be enhanced by spillover, due to the strong affinity of H atom on COFs surface, the high potential energy is required in the dissociation process. Therefore, it is necessary to discover other catalysts to reduce the migration and desorption barriers of H atoms.
Actually, the metal catalyst plays a very significant role in the spillover process; for example, to a great extent, it may determine whether the dissociation of H2 molecules and migration of H atom toward the substrate can occur. Therefore, it is of great importance to study the effects of different metal catalysts on the entire spillover process in COFs, however, to our knowledge, such effects have not been reported to date. Currently, the most widely used metal catalysts in the spillover process are Ni, Pd, and Pt nanoparticles. Therefore, this paper aims to investigate the effects of different metal catalysts (Pt, Pd, and Ni) on hydrogen spillover mechanism and evaluate the feasibility of hydrogen spillover on the COF-105 via DFT calculations. We choose COF-105 as the model system because it is a typical three-dimensional (3D) COF with the highest hy drogen storage capacity among all the COFs, i.e., 21 wt% and 4.5 wt% at 77 K and room temperature with the pressure of 100 bar, respectively.[10] Herein, we mainly study the depositions of metal catalysts M4 (Pt4, Pd4, and Ni4) on COF-105, adsorptions and dissociations of H2 molecules onto metal catalysts, H migrations from metal catalysts toward COF-105, and surface diffusion of H atom. The effects of Pt4, Pd4, and Ni4 catalysts on each process are also compared with each other in detail.
2,3,6,7,10,11-hexahydroxy triphenylene (HHTP) and tetra (4-dihydroxyboryl-phenyl) silane (TBPS) clusters were utilized to represent the building units constructing the COF-105 in order to reduce the computational cost, which were cut from the optimized periodic structure and terminated the breaking bonds with H atoms as shown in Figs.
All DFT calculations were carried out by using the Vienna ab initio simulation package (VASP)[24,25] with the exchange–correlation function of Perdew–Burke–Ernzerhof (PBE).[26] The core and valence electrons were described by projector augmented waves (PAW)[27] and a plane-wave basis set, respectively. An energy cutoff of 500 eV was used to make sure that the total energy is well converged. To avoid interactions between adjacent periodic images, the cluster models investigated in this paper were placed within a 25 Å × 25 Å × 25 Å box. The structural optimization was performed with the energy and force converging to less than 10−4 eV and 0.03 eV/Å, respectively. A 2×2×2 Monkhorst–Pack mesh was used to represent the Brillouin zone.
The Climbing Image Nudged Elastic Band (CI-NEB) method [28] was used to determine the minimum energy pathways of migration and diffusion of H atoms. All of the initial and final states were from optimized configurations. Six images between initial and final states were chosen to obtain smooth potential energy curves.
First, we investigate how metal catalyst clusters M4 (Pt4, Pd4, and Ni4) are deposited on COF-105. Three deposition configurations are considered, which are point, edge and surface of M4 contacting with HHTP and TBPS clusters, respectively. As a representative, the point, edge and surface deposition configurations of Pt4 on HHTP and TBPS clusters are shown in Fig. (color online) Schematics of optimized deposition configurations of Pt4 on HHTP and TBPS clusters, where cfg1, cfg2, and cfg3 represent point, edge, and surface configurations, respectively.
The binding energies of point, edge and surface configurations of M4 cluster on the substrate are listed in Table
The binding energies of M4 on HHTP and TBPS are 1.65 eV–1.88 eV and 2.26 eV–2.57 eV for Pt4, 1.26 eV– 1.79 eV and 1.92 eV–2.34 eV for Ni4, 0.53 eV–1.29 eV and 1.69 eV–2.04 eV for Pd4, respectively, which indicates that among the interactions of Pt4, Pd4, and Ni4 with the substrate, Pt4 is the strongest, Ni4 is the middle and Pd4 is the weakest. This order is consistent with that of Sigal et al.[29] They used the DFT method to investigate the adsorption energies of hy-drogen and oxygen on graphene decorated with a wide set of metals (Li, Na, K, Al, Ti, V, Ni, Cu, Pd, Pt). In conclusion, M4 cluster is preferentially deposited in the vicinity of TBPS cluster and tend to be in the form of surface contact; Pt4 and Pd4 have the strongest and weakest interactions with the COF-105 surface, respectively.
In hydrogen spillover process, H2 molecule first dissociates into H atoms around the metal catalyst through the interaction with metal catalyst, then the dissociated H atoms migrate toward the substrate and diffuse in the whole substrate surface. Thus the key conditions for hydrogen spillover to take place are as follows: 1) the H2 molecules adsorbed onto metal catalyst can dissociate into H atoms; 2) the dissociated H atoms can overcome low energy barrier to migrate toward the COF-105 substrate; 3) H atoms can further diffuse in the COF-105 surface.
In this section, we consider the adsorption and dissociation of H2 molecules on the M4 catalyst. The saturated hydrogen adsorption configuration on the M4 is the initial state of hydrogen spillover. Thus we successively place H2 molecules in the vicinity of M4 for M4@substrate. Without the contact between M4 and COF-105 surface after the optimization, up to 16 and 14 H atoms can be chemisorbed on Pt4 for Pt4@HHTP and Pt4@TBPS, respectively, which indicates that the dissociation process of H2 molecules is automatically to meet the first condition of hydrogen spillover. This result is consistent with that of Ref. [30]. The full H saturated adsorption configurations of Pt4 for Pt4@HHTP and Pt4@TBPS are shown in Fig.
The average chemisorption energy of H2 molecule on M4 cluster is defined as
The average chemisorption energy ∆EChem values of H2 on M4 clusters at full saturation with different configurations are obtained according to Eq. (
The effects of Pt4, Pd4, and Ni4 on ∆EChem are also compared. It can be seen that from Table
The migration process of H atoms from M4 toward the substrate surface is well studied by using the CI-NEB method. The adsorption locations of H atoms on M4 clusters can be divided into the top site (adsorption at the top of metal atom in the form of metal–H bond) and the bridge site (adsorption onto the metal–metal bond in the form of metal–H–metal bond). The results show that only the H atoms on the bridge site can migrate to the COF-105 surface. However, for the H atom on the top site, even if we dragged it onto the C atom in the substrate, it would rebind to the top site of M4 after optimization. This is mainly because the interactions of metal atoms with the top H atoms are stronger than those with the bridged H atoms, thus the H atoms onto the top site do not migrate easily. Wu et al.[32] also found a similar phenomenon when they investigated the migration of H atom from Pt4 clusters toward B-doped graphene surface. It is noted that, when the H atom located on the bridge site (labeled as 1) migrates to the substrate surface, the other H atom located on the top site (labeled as 2) will migrate to this bridge site, resulting in the migration proceeding continuously and satisfying the second condition of hydrogen spillover. Especially for the H atoms onto the Ni4 cluster, the activated H2 molecule will be further dissociated into H atoms (labeled as 3 and 4) (see the inserts in Figs.
As surface contact configuration is most stable for M4 deposited on the substrate as discussed in Subsection 3.1, we take the cfg3 configuration for example, and use the CI-NEB method to study the minimum energy pathways of the migration of H atom from M4 toward the HHTP and TBPS surface, as shown in Fig.
The diffusion processes of H atoms on M4 supported and pure COF-105 surface are investigated, respectively. Their minimum diffusion pathways are shown in Fig.
The effects of Pt4, Ni4, and Pd4 catalysts on the diffusion process are also compared. As shown in Fig.
Combining the above discussion for the whole spillover process, we conclude the relationship between the metal catalyst and each step of spillover process. The binding energy of M4 on the substrate
The hydrogen spillover mechanism on metal-supported COF-105 is investigated by DFT calculations, in which included mainly are various deposited configurations of metal catalyst on the COF-105, adsorption and dissociation of H2 on the metal catalyst of M4@COF-105 migration of H atoms from M4 toward the COF-105 and diffusion of H atoms on the COF-105 surface. The results are shown as follows.
The larger the number of metal atoms interacting with the substrate, the greater the binding energy of metal deposited on it will be. Metal catalyst prefers to deposit on TBPS cluster with surface contact configuration. Only H atoms located at the bridge site may migrate to the substrate surface, however, the migration process is an endothermic reaction and less stable. The introduction of metal catalyst can greatly reduce the diffusion energy barrier of H atoms on substrate surface, which makes it easier for the H atoms to diffuse on the substrate.
The effects of metal catalysts M4 (Pt4, Pd4, and Ni4) on the hydrogen spillover process have also been systematically compared. It is found that the order of binding energy of M4 deposited on the substrate is Pt4 >Ni4 >Pd4. The migration and diffusion barriers of H atom are lower for Pt4 catalyst than those for Ni4 and Pd4 clusters, which indicates that Pt4 is the most favorable for hydrogen spillover on COF-105. However, the H migration process from Pt4 toward the substrate is still endothermic. Therefore, it is thermodynamically unfavorable that the H atoms directly migrate from metal catalyst toward the substrate. More efforts are needed to further improve the thermodynamical feasibility of hydrogen spillover, which is what we are going to do in the next step.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] |